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We analyse the possibilities for quantum state engineering offered by a model for Kerr-type non- 
linearity enhanced by electromagnetically induced transparency (EIT) , which was recently proposed 
by Petrosyan and Kurizki [Phys. Rev. A 65, 33833 (2002)]. We go beyond the semiclassical 
treatment and derive a quantum version of the model with both a full Hamiltonian approach and 
an analysis in terms of dressed states. The preparation of an entangled coherent state via a cross- 
phase modulation effect is demonstrated. We briefly show that the violation of locality for such 
an entangled coherent state is robust against low detection efficiency. Finally, we investigate the 
possibility of a bi-chromatic photon blockade realized via the interaction of a low density beam of 
atoms with a bi-modal electromagnetic cavity which is externally driven. We show the effectiveness 
of the blockade effect even when more than a single atom is inside the cavity. The possibility to 
control two different cavity modes allows some insights into the generation of an entangled state of 
cavity modes. 

PACS numbers: 03.67.-a, 42.50.Dv, 42.50.Gy, 42.65.-k 



I. INTRODUCTION 

The reliable preparation of non-classical states of light 
such as a travelling- wave entangled coherent state P, 13 
and Schrodinger cat state , or the control of the popu- 
lation of an individual field mode (of an electromagnetic 
cavity, for example) are recognized to be important tasks 
in Quantum Information Processing (QIP) jj]. Entan- 
gled coherent states and Schrodinger cats, for example, 
revealed useful for QIP with coherent states On the 
other hand, photon blockade appears as a striking man- 
ifestation of control on a system at the quantum level 
and opens a way to novel schemes for quantum state en- 
gineering 5] . 

In this paper, we investigate the possibilities of quan- 
tum engineering using non-linear processes realized by 
electromagnetically induced transparency (EIT) The 
EIT regime seems to be able to overcome one of the ma- 
jor bottlenecks in a non-linear process: the low efficiency 
accompanied by high absorption rate of a conventional 
Kerr medium, that makes the production of a travelling- 
wave cat state, for example, far from realization 0. 

On the other hand, it has been proved that the atomic 
medium in EIT regime shows a measured x^^^ parame- 
ter up to six orders of magnitude larger than usual @. 
This suggest the use of the enormous non-linearity to 
get a reliable cross phase modulation effect between two 
travelling fields of light even for the very low photon- 
number case 0, 01 . Usually, the approach to such a 
process is semiclassical. But, for the purposes of QIP, a 
full quantum treatment is relevant & illj . Such an anal- 
ysis has been performed in ref. where the idea of 
a double-EIT regime is introduced in order to optimize 



the non-linear interaction between two electromagnetic 
fields. The proposed model has been modified, in (ISj . to 
get an easier experimental realization of the process. In 
this latter work, however, the analysis is again semiclas- 
sical. We give the full quantum mechanical description 
of 13] by means of an Hamiltonian approach 14] • 

A promising candidate for the embodiment of the 
atomic model we discuss is a Pr'^^ doped Y2Si05 crys- 
tal (Pr:YSO): it has an energy-level scheme appropriate 
for our purposes and it has been used for the demonstra- 
tion of EIT 15] and giant non-linearity in solid state de- 
vices [l6|,[lj|. Using typical values for Pr:YSO, we find a 
giant non-linearity even at the quantum level. We derive 
the equations of motion for the involved quantum fields 
from an effective interaction Hamiltonian. With these re- 
sults, starting from two independent coherent states, we 
prepare an entangled coherent state and a Schrodinger 
cat state. 

As a second example of the applicability of these 
results, we treat a cavity-quantum-electrodynamics 
(CQED) system. We concentrate on a photon-blockade 
effect realized combining the obtained large non-linearity 
and the features of isolation from the environment and 
manipulability characteristic of CQED. As we see, the 
interaction of the atomic model we depict with the elec- 
tromagnetic field of a cavity results in a coupled system 
that exhibits a non- linear eigenspectrum This is ex- 
poitable to control the number of excitations that are fed 
into the cavity from an outside radiation. Furthermore, 
the specific model for a double-EIT allows to treat the 
interaction with a bi-chromatic cavity field and to show 
how to manipulate it to settle entanglement. 

The paper is organized as follows: in Section ^ we 
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sketch the Hamiltonian method we have chosen; it is im- 
mediatly applied, in Section [llll to the atomic model for 
double-EIT We then derive the equations of mo- 

tion for the quantized fields and give the order of mag- 
nitude of the achieved rate of non-linearity. Section IIVI 
is devoded to an alternative approach to the double-EIT 
problem: we choose a dressed state picture to re-derive 
the polarizabilities of the medium. The effective interac- 
tion Hamiltonian derived in Section llTTl is used to show, in 
Section[V] how a tensorial product of two coherent states 
evolves toward an entangled coherent state. A technique 
to project one of the modes onto a Schrodinger cat state is 
analyzed. We describe a scheme for the characterization 
of a generated cat state . In Section IVII we give the 
outlines of a bi-chromatic photon-blockade realized, in an 
electromagnetic cavity, by the high rate of non-linearity 
inherent in the chosen model for double-EIT. 

II. THE HAMILTONIAN METHOD 

Usually, the interaction of electromagnetic fields with 
an atomic medium is mathematically described by means 
of the Maxwell-Bloch equations, which are a set of cou- 
pled differential equations that connect the dynamics of 
the fields to that of the atomic degrees of freedom. These 
latter evolve according to the Von Neumann equation 
ih ^'^gP = [H,p{t)], with p{t) the atomic density opera- 
tor and H = Hatom + H field + Hinteraction the COmplctC 

Hamiltonian model for the problem. Whenever condi- 
tions of adiabaticity and low intensities of the fields are 
valid, a steady state solution can be obtained. Inserting 
it into the Maxwell equations gives the evolution of the 
fields. 

In an EIT problem, however, the number of interacting 
fields as well as the involved atomic levels is usually large. 
This makes the analytical solution of the Maxwell-Bloch 
equations a challenging task. From classical considera- 
tions, it is possible to see that the polarization Pj of a 
medium can be expressed as: 

Nd.j I dH'\ , . ,x 

)+c.c., (1) 

where Q.j is the Rabi frequency of the j-th field (of fre- 
quency ojj), dj is the dipole matrix element of the cor- 
responding transition, N is the density of the atomic 
medium and H' is the single-particle interaction Hamil- 
tonian. Here we assume that the atoms in the medium 
are equally coupled to the different fields. When Eq. ^ 
is introduced into the Maxwell-Bloch equations and we 
use the slowly varying envelope approximation (SVEA), 
we get: 

f d I Nd^Mj I dH'\ ^ 

The fields responsible for the coupling of the initially 
prepared atomic state to other levels are taken, here, to 



be weak (weak coupling limit). This gives a small proba- 
bility of transition toward states different from the initial 
one. The initial state becomes, thus, a kind of stationary 
state, whose evolution will be adiabatically followed by 
the ensemble. In Eq. jSJ, H' can thus be replaced by A, 
the energy eigenvalue of the initially prepared state. 

The quantization of the interacting fields can be done 
giving operatorial nature to the field variables in the ef- 
fective Hamiltonian that A represents. Bosonic commuta- 
tion rules to the fields creation and annihilation operators 
are thus imposed |T^. 

III. CROSS PHASE MODULATION VIA A 
DOUBLE EIT EFFECT 

Here we describe the model for double-EIT proposed 




FIG. 1: The atomic model for double-EIT. Fields Ea and 
Efj have frequency uJa and oJb, respectively. The fields with 
frequencies LUd\, (^d2 are classical (their intensities are much 
greater than those of Ea and Eb). Detunings Au 7^ A_l and 
jA| = \Au — Ai|. The decay rates of the excited states are 
taken equal to 7, for the sake of simplicity. 

It involves a ground metastable triplet {|1) , |2) , |3)} 
and an excited one {|4) , |5) , |6)}. By means of an exter- 
nal magnetic field the ground states are splitted by 
while Ajj (7^ Al) is the splitting for the excited triplet. 
Two weak fields, Ea and Et, drive resonantly |2) |4) 
and 1 2) ^ |6) respectively; |1) |5) and |3) ^ |5) are in 
the dispersive regime (detuning |A| = \Au — Al|) while 
transition |2) ^ |5) is assumed to be forbidden. The 
couplings |1) ^ |4) and |3) ^ |6) are realized by two 
classical fields of different frequencies but equal Rabi fre- 
quencies. Two distinct subsystems are easily singled out: 
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states , |4) , |2} , |5)} constitute the four-level N sys- 
tem proposed in Q to show giant Kerr non-linearity. EIT 
is realized for field Ea while Ef, acts as a perturbation 
that induces a.c. Stark shift on |1), which determines 
a large non-linear effect An analogous description 

can be done, interchanging Eb with Ea, for the subset 
{|3) , |6) , |2) , |5)}. The two subsystems are connected by 
the non-resonant couplings to 1 5) . In the interaction pic- 
ture, the Hamiltonian reads: 



H' = h{{A\5) (5| + 

+r!f,|5)(i| + r!fc 



QdH) (l\+nd\6) (3\+fla\4:) (2| 
6) {2\ + Qae-^'^' \5) (3| +C.C.}. 



(3) 

We then introduce the decay rates of the excited states 
7 and phenomenologically change the signs in front of 
each Rabi frequency in Eq. ijjjl to match the expression 
in ref. 13] . The atoms in the ensemble can be prepared 
in state |2) by means of optical pumping and, if the weak 
field limit is assumed {\^d\ ^ l^^a,b|), the atomic system 
remains in |2) all along the interaction time. 

To discard the Doppler broadening, an atomic gas can 
be mantained at a low temperature. Alternatively, we 
can take a solid state device as the Pr^+ doped ¥28105 
crystal (Pr:YSO) [23| to embody the Hamiltonian model. 
This latter choice is motivated by the similarity between 
the energy-level scheme described here and that of the 
transition ^^^4 — > in Pr:YSO. In ultraslow 

group velocity (~ 45 m/sec) in Pr:YSO has been re- 
ported. 

Following the recipe outlined in Section ^1 we solve 
the secular equation for H' and seek for the eigenvalue 
of the initial state |2). In the weak field limit we get: 



2h\na\^\nb\^ 



(4) 



The derivative of Eq. Q with respect to fi* (fi^ ) gives 
the polarization of the medium at frequency tOa {(^b)- The 
equation of motion for f2a can be, finally, cast into: 



d_ 

dx 



ld_ 

cdi 



2zjV(7o7|»fcr 



rfla = icXa^a, (5) 



with aa the atomic polarizability at frequency LUa |l3j| . 
Here, ctq = ■jeoch'y resonant absorption cross sec- 

tion. The equation for Eb can be analogously derived. 

It has been shown in ^3 that the interaction between 
two fields in a medium exhibiting EIT is critically lim- 
ited by the time that the faster of them spends inside 
the medium itself. This velocity mismatch strongly af- 
fects the efficiency of any non-linear process we want to 
realize in the usual N configuration. In (12i |. the induc- 
tion of an EIT regime for both the fields (double-ElT) 
is suggested to bypass the problem. Strongly reducing 
the group velocities of the beams (by means of EIT), 
the interaction time could be maximized, optimizing the 
efficiency of the non-linear process. If we compute the 



group velocities for fields Ea and Eb in the model de- 



scribed above, we find ■ 



\ndr/Nao^ < c. This is 



the signature of the double-EIT established in the atomic 
ensemble. 

We now demonstrate how this effect is useful for a 
strong effective non-linear effect on the evolution of Ea 
and Eb in the full quantum domain. We give operato- 
rial nature to the Rabi frequencies in Eq. I^J introduc- 
ing the positive and negative frequency components of 
the corresponding operators. For example: ria{z,t) — 

^^24Efcl/5S^«fe(*)e"*^'^''""'"^^'^' (and analogous for 

Clb). This expression well describes a pulse in the nar- 
row bandwidth approximation: fc is a label for the wave- 
lengths in the packet, w^'"' is the central (carrier) fre- 
quency of the pulse and Vq is the quantization volume. 
The narrow bandwidth approximation takes Slo <C lu^'^^ , 
where Slo is the bandwidth of the pulse. These operators 
satisfy the commutation rules [Cti, fij] cx Sijl — a, b), 

with 6ij the Kronecker symbol and 1 the identity opera- 
tor. Multiplying by the atomic density N and integrat- 
ing over the interaction volume V, we have the effective 
Hamiltonian: 



Heff = 



2hN 



nlCianldb 



{ii-A)jy \nd\ 



dV. 



(6) 



We can now write the Heisenberg equations for Ea and 
Eb- The evolution of the probe fields is, then, given by: 



EaAL,t)^ Ea.b{0,t')e 



(7) 



with t' = t — L/v^^°'^P, X a non-linearity rate (obtained 
collecting all the non operatorial quantities in Eq. ©) 
and L the interaction length. The above equation shows 
explicitly an effect of cross phase modulation on the 
quantum fields. For a cw laser beam, we can recast 
Eq. © as: 



TtCW 

^eff 



with X ' 



hxa^ dwh 

NuJaUJb M24|^ M26|^ 

2tfel{i-i-A)\nd\^V 



(8) 



where we have assumed V = Vq and 5R{-} is to take the 
real part. We can write and solve the Heisenberg equa- 
tions for a and 6, representative of the dynamics of Ea 
and Eb respectively, getting: 



-iXtfii, 



a(0) 



tit) 



-iXtfla 



6(0), (9) 



where fij is the photon number operator for field j — a, h. 

The evolution of one field rises up depending, in the 
most explicit way, on the intensity of the other one. To 
estimate the order of magnitude of x we use values for 
the parameters in (jSJ typical of the ^H^ transi- 
tion in Pr:YSO. We take a wavelength of ^ 600 nm [l^ . 
L ~ 1 mm and a cross section for the beams of 100 ^m 
(FuU-Width-at-Half-Maximum) . The decay rate 7 can be 
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taken between 10 and 100 kHz (sample's temperature of 5 
K); A ~ 1 MHz and \ flci\ ^ 1 MHz are reasonable values 
and allow us to consider 7 ^ A, \ild\, which gives a negli- 
gible rate of two-photon-absorption. Indeed, two-photon 
absorption appears because of the two-photon resonant 
driving at the basis of EIT. The absorption is propor- 
tional to the imaginary part of the polarizability of the 
medium and it can be minimized if A ^ 7. The elec- 
tric dipole matrix elements for the system in exam are 
typically - lO'^^ Cm and N ~ 10^^ cm'^. With these 
values in Eq. |(SJ) and for interaction times T ~ /isec, a 
cross-phase shift xT = tt is achieved, even with beam 
intensities of just a few of photons. 

An unwanted effect comes from the, here neglected, 
couplings of Ea, Eb with the atoms via the transitions 
|1) ^ |5) and |3) ^ |5) respectively. These spuri- 
ous couplings lead to the self-phase modulation, an ef- 
fect for which a field evolves independently from the 
other one. The polarizability, in this case, scales as 
ctj'^'''^ oc ftj / {Ajj + Al). Because this is smaller than the 
effect of cross phase modulation and independent from it, 
we can safely neglect it when dealing with cross effects. 



IV. ATOMIC POLARIZABILITIES IN A 
DRESSED STATES PICTURE 

In what follows, we change the point of view and we 
treat the problem of double-EIT using a Hamiltonian 
model in which the fields are quantized ab initio. We 
refer again to Eq. ^ but we now make the subtitutions 
na |4) (2| ^ gaa |4) (2|, fti, |6) (2| ^ gtb |6) (2|, where ga,b 
are related to each coupling. Clearly, we have to add, to 
Eq. the term hiojafia + tUbfib) that is the energy of 
the quantized fields. We, thus, can write: 



H' ^h{nd\4){i\-i^d\6) (3| 

+gaa\4) (2|-5b&|6) (2| +c.c. 



(10) 



where the couplings |1) ^ |5) and |3) <-> |5) have not 
been introduced, yet. We justify this choice later in this 
Section. It is possible to show that the states coupled by 
the Hamiltonian H10|l are the elements of the following 
set: 

{\l,na - l,nb) , \2,na,nb) , \3,na,nb - 1) , 

|4, Ua - 1, Ub) , |6, Ua, Ub - 1)}. 

The matrix representation for iJ', in the Hilbert space 
restricted to these states, is: 



with the Rabi frequencies taken as real. The diagonal- 
ization of (|11|) leads to the eigenvalues: 



E,^0 



E{ = -hn E+ = nn 



(12) 



e: 



-Mlrl 



Mid 



with n = + gl{na) +gl[nb). The eigen- 

states can easily be found and we call them 
{|0«„„.) , |1- „J , |1+ „,) , \K.n,) , |2+ „,)}. This 
dressed eigensystem shows that the model can be 
mapped into a five-level ladder model. This map is 
just a rotation of the bare basis into the dressed one: 
I dressed) = Pjbare), where P is the matrix that 
realizes the rotation; | dressed) and |bare) are two 
vectors of dressed and bare states, respectively. 

On the dressed states, Eq. (|10|l is diagonal and takes 
the form: 



^'dressed — ^ ^ ^ 
1=1.2 j=± 



(13) 



We introduce the couplings to level |5) writing the spin- 
flip operators |5) (1| and |5) (3| in terms of dressed states. 
To do it, we use a closure relation so that, for example: 



. ^d9a^/K 



i; 



^dga^a /-,+ 



9b ^ 



with 6' = 



■9b<. 



V2S' 

(14) 

Analogous expressions can be 
derived for the other field-atom operators. Furthermore, 
the term A |5) (5| has to be added to H10|l . This way to 
introduce the bare atomic level |5) reminds a perturbative 
approach justified because of the dispersive nature of the 
couplings to this level. In the weak field limit and for a 
sufficiently large detuning A, the transition probability 
to 1 5) remains small. 

To shorten the notation, in the following we take 
9a\pn^ = f^a, 9b\/nb = ^b- Applying the Hamiltonian 



operator to a generic state vector, decomposed as: 











^d 













9ay/n^ 


-9bVnb 
















^d 


9ay/nZ 











Vo 


-9b^/nb 


~^d 





/ 



(11) 



1^) = Ao |0„„„J + |1- + A+ |1+ „J 
+ A^ |2-„„J + A+ |2+ „J + ^5 \5nanb) , 



(15) 



gives a Schrodinger equation that is equivalent to the 
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following set of differential equations: 



idtAa = 


\ f, ^\ f.i, 

n 


A, 


idtA^ = 


EiA^ 




idtA^ = 


E^A^ 




idtA+ = 


EtAt 




idtA+ = 


EtAt 


A, 



(16) 



idtA, = AA5 + ^ {A+ - A^- 



We have neglected terms oscillating at frequency 2A + 
E^2 because they average to zero when the time integrals 
are carried out. Here, again, we adopt SVEA and we 
assume as initial state |2). The contribution of this bare 
state to the linear combinations that define the dressed 
ones is relevant only for |0„^„J, that has the form 



|0„„„J = - ^ |1, - 1, Ub) + ^ |2, Ua, Ub) 
- ^ |3, na,nb - 1) , 



(17) 



while it is of order (5/51 ^ 1 or even null for all the other 
dressed states. Thus, for il^ — fl, we can take ^0 = 1- 
Moreover, it is easy to show that A'^{t) = A^{t) and 
A-{t) = A^{t) which leads to A^it) ~ jiaf^{,/(^^A). We 
note that |0„^„(,) is a dark state since it is composed of 
the ground states only and does not contain the decay- 
ing states 1 4) and |6). All the other eigenstates have a 
contribution from both |4) and |6) and are bright ones. 
We can derive all the other probability amplitudes: 



A^ 



^g^b^ri A+ = A^ = ^"-^b 



(18) 

If we want the atomic polarizability at frequency uja, 
as in Eq. lO, we have to come back to the bare state 
description. To find the probability amplitudes for each 
bare atomic state, we write = X^^i A I*) equate 
it to its expression in terms of dressed states. We look 
for the coefficients 02, Ps, (34, P5, since the polarizability 
aa is proportional to {^2^4 + ^3^5)- Clearly f]^ = A5 
and for the others, we have: 



P2 



^d 

n 



Pi = (4|V) = 



03 



im = 



n 



(19) 



so that aa — ^'^n^^'' ' with (Tq as defined in Section IlIII 

d 

This result matches what has been obtained by the 
Hamiltonian approach once we replace the field variables 
with the corresponding operators. The dressed states ap- 
proach described here and the Hamiltonian one lead to 
consistent results. Since the former relies on an ab initio 
quantum level, an undeniable reliability is given to the 
latter method. 



V. SCHRODINGER CAT STATES 
GENERATION 

We apply the results obtained in the full quantized 
picture of the cross phase modulation via double-EIT to 
produce non-classical states of field mode. The evolution 
shown in Eq. can be derived from the action of the 
unitary operator U{(p{t)) — g^-^vWriant ^ .^^j^^j^ ^^(^^^ _ 
PH . If the initial state of the field modes E^, Eb is 
the tensorial product of two coherent states \'4'{Q))ab ~ 
|a)^® 17)^, its evolution by means of U{(p), for (p{T) = n, 
is given by 0| 

m^/x))ab « l«)a {l7) + \-l)}b + |-«)a ill) " hl)}b ■ 

(20) 

This is a particular expression for an entangled co- 
herent state: unitarily acting on the subsystem b we 
can transform it into the more familiar form \a)^ \"f)f^ + 
\—c()a \^l)b- ^1- the superpositions of coherent 
states 1 7) J and |— 7)^ are Schrodinger cat states: 



l7>, 



-7)6 a 



7 



2j 



l7)b - |-7)b « 



2j + l 



|2j>b, 



(21) 



7 



|2j + i; 



b ■ 



These are sometimes called even and odd coherent states. 
We have shown how the described non-linear interac- 
tion, for the particular case of initially prepared coherent 
states, results in an entangled state. 

To prove the entanglement, we have to show the cor- 
relation of the modes a and b as we unitarily transform, 
gradually, from [7)^ to |7)j, ± |— 7)},. However, this in- 
volves another non-linear interaction. Wc thus discuss 
an indirect procedure. In details, if we can discern where 
field a is, the state of field b is projected onto one of 
|7)j ± |— 7)},. To determine the state of field a we use 
a 50 : 50 beam splitter (BS) and two photodetectors, 
as sketched in Fig. |21 After passing through a BS, two 
coherent input fields \a) , become: 



a + l3 



V2 



-a + P 



V2 



(22) 



where 5a 



^ (a^ c—dc^ ) 



is the BS operator, with a (a) 



and c (c) its input (output) modes. For P = a we have 
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the following read-out: if the input mode a is prepared 
in \a)^, then Detector 1 will reveal some photons, while 

Detector 2 will not. In this case, the field mode b will 
be projected in the even coherent state. In the oppo- 
site event, mode b will be in the odd coherent state. Of 
course, there is a possibility to have both the detectors 
not to click: we do not know in which state mode a is 
and we have to repeat the procedure untill one of the 
detectors clicks. 



field mode a 



Detector 1 
field mode a 




Detector 2 



FIG. 2: Scheme to infer the state of the field mode a. It is 
shown the symbol used for a photodetector and that for the 
50:50 beam splitter (BS). Detector 1 (2) clicks just if the state 
of mode a is \a)^ {\~'^)a)- With this scheme, we can generate 
an even or an odd coherent state of mode 6. 

A possible way to detect the quantum nature of the 
state generated by the described scheme is schematized 
as follows: to indirectly infer the coherences in an even or 
odd coherent state we generate a new entangled coherent 
state, mixing field mode b (which is supposed to be in one 
of the states in Eq. (|21|l 'l and the vacuum of an auxiliary 
mode in a 50 : 50 BS. Calling b,c the output modes of 
the BS, the resulting joint state of radiation is 



4 = -^ 



7 -7 \ / 7 -7 



"7 J_\/zl _L 
^/2' V2/\^/2'^/2 



-7 7 
V2' V2 



7 -7 
^/2' ^/2 



7 -7 
V2' V2 



^/2'^/2 



be 

(23) 

where A is a normalization constant. The parameter c 
takes account of the purity of the generated Schrodinger 
cat state. If it is c = 1 (c = — 1), the field mode b was in 
an even (odd) coherent state, while for c = 0, it was a sta- 
tistical mixture that has been produced by the non-linear 
interaction. To discern between the possible values for c, 
we use the criterion for inseparability proposed in 
and we evaluate the function S = ((Au)^) -t- ^(Aw)^), 
with (Am)^ and (A-0)^ the variances of tt — + Xc, 
V = —pi + Pc and {xj,pj} the phase-space quadrature 
operators for mode j = b,c . 



According to the sufficient condition for inseparabil- 
ity in 0, if 5 < 2, the state of b and c is inseparable. 
To experimentally evaluate 5' we need the single quadra- 
ture variance {(^{Axa.b)"^) for example) and correlations 
as (xaXb) or (paPb)- AH these quantities can be deter- 
mined via two homodyne detectors, one for each field 
mode [23, ■ A plot of 5 as a function of the ampli- 
tude 7 in the case of an even coherent state is given in 
Fig. 121 (a). The bound 5* = 2 is violated, for c = 1 just 
until 7 ~ 2 (we have S = 1.995, for 7 = 2), revealing 
the entanglement of the generated state. Our criterion 
here is only a sufficient condition for entanglement which 
works fine for a small number of photons. As 7 grows, 
this sufficient condition does not give information on en- 
tanglement. 




FIG. 3: (a) Plotof S = {{Auf) + {{Avf) as a function of the 
amplitude 7. The effect of different values of the parameter c 
in the density matrix p'^ . is studied: the dot-dashed curve is 
for c = 0, corresponding to the case of a statistical mixture. 
The dashed curve is for c = 0.5 while the solid curve represents 
a perfectly generated even coherent state, (b) Behaviour of 
the separability function when imperfection in the homodyne 
detection is considered. Here, the dot-dashed curve is for 
detection efficiency rj = 0, the dashed one is for rj = 0.4 and 
the solid curve is for r; = 0.8. For lower jy, the minimum values 
of S feet shift toward higher 7 values: this is because the 
lower is the efficiency of the homodyne detectors, the more 
the input state resembles a Gaussian state. 



It is possible to introduce the homodyne detectors 
losses modeling an inefficient homodyne detector with a 
beam splitter BS,; (transmittivity 77) followed by a perfect 
homodyne detector 0]. Each beam splitter BS,, mixes a 
mode of the signal to measure with a vacuum state and 
transmits the signal with probability 77; the amount of 
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reflected input field is a measure of the losses. The quan- 
tum efficiency of the detectors is, thus, rj. The calcula- 
tion of the total variance for the quadratures of modes 
6, c when the detectors have an equal quantum efficiency 
77 leads to what is shown in Fig. (b). The separabil- 
ity function keeps its functional features even in the case 
of lossy detection and some similarities with the case of 
perfect detection are evident, showing the robustness of 
the scheme. 



VI. BI-CHROMATIC PHOTON BLOCKADE 

In this Section we give the outlines for a possible quan- 
tum control of light offered by the giant non-linearity, in 
the quantum regime, of the model for doublc-EIT. In par- 
ticular, we analize the interaction of a low density group 
of atoms, with the energy scheme sketched in Fig.^ with 
two modes of an optical cavity. Each cavity mode takes 
the place of one of the probes in Section Hill and the cav- 
ity itself is driven by two weak external beams, each one 
on-resonance with a relevant cavity mode. The driving 
fields at frequencies LUdi , '-^d2 are then shined on the atoms 
to obtain the double-EIT regime. The physical system is 
schematically shown in Fig. ^ 




FIG. 4: Set-up for a bi-chromatic photon blockade via large 
non-linearity. The optical cavity is crossed by a low density 
beam of atoms, each one having the six-level energy scheme 
suitable for a double-EIT regime. This condition is estab- 
lished by two classical driving fields (not shown in the pic- 
ture) and by two cavity field modes, driven on-resonance by 
two external fields, £2""* and 

For the moment, we treat the case in which the density 
of the atomic beam is so low that the cavity is crossed by 
a single atom, each time. In this condition, the Hamil- 
tonian of the system atom+ cavity + external fields is H — 
H' + h jwao^a + Wf,6t5| + nSpump {(a^ +a) + {b^ + 1)^ 

where H' is the Hamiltonian in Eq. H1(J|I and the third 
term takes account for the coupling of the cavity with 
the external fields. The parameter Epump is taken the 
same for both E'^^* and i?^^*. The same arguments used 
in Section 1^ lead to neglect, for the moment, the cou- 
plings to the atomic level |5) while, in general, terms for 
the damping of the cavity modes have to be added. We 
will consider these two points later. As before, the initial 
state for the atoms is |2). 



We examine what happens if E^^* shines the cavity, 
exciting the corresponding cavity mode. In this case, the 
six-level atomic model reduces to a simpler N configura- 
tion, where the perturbation to the otherwise perfect EIT 
regime for Ea is due to the |1) <^ |5) couphng by Ea itself. 
The resulting non- linear self-phase modulation effect, de- 
scribed in Section IIIII can not be neglected because the 
cross phase analog is not active. A mono-chromatic pho- 
ton blockade effect results, as studied in Ref. |5l l23.l2q . 
Once a photon in E'^* leaks into the cavity and feeds 
Ea^ no other external photons are allowed to enter (the 
photons are blocked). This is because the self-phase ef- 
fect gives non-linear features to the atom-|-cavity mode 
system: in a dressed-state picture, a second external pho- 
ton is resonant with no one of the transitions that lead 
from a singly excited dressed state to a doubly excited 
one. When the mono-chromatic photon blockade is ac- 
tive, the dressed system atom-|-cavity mode is trapped 
between the bare state |2,0)^jq„j^ and the dark state 

\Da) = N {nd\2,l) - ga\l,Q)}atol^E^ Mh a normal- 
ization factor M) and behaves as an effective two-level 
system. The arguments above can be reformulated when 
j^jext fgg(jg lY^Q cavity. 

We now treat the situation in which the cavity has 
been already fed by one photon (for example from i^^^*, 
so that Ea is excited) and we consider the effects of the 
interaction with a photon of different colour (i.e. with a 
photon from E^^*). The initial state is |2, 0, 0)atom Ea Sb' 
that is the state without excitation. The absorption of 
the first photon by the cavity takes the system to the 
dressed state \Da). To see if a photon blockade with re- 
spect to incoming photons of frequency cob is possible, we 
diagonalize the interaction Hamiltonian in the subspace 
spanned by states as \atom) (g) \Ea) (8) \Ei,) that have two 
excitations. 

The interaction Hamiltonian, taking explicitly into 
account the couplings to |5), is closed within 
{|2, 1, 1) , |4, 0, 1) , |1, 0, 1) , |6, 1, 0) , |3, 1, 0) , |5, 0, 0)}. We 
take A fid and Qa ^ gb- As a function of the Rabi 
frequency gb, a typical plot of the eigenenergies is shown 
in Fig. (a). A photon blockade effect is evident: to 
feed the cavity, a photon of frequency lo^ should find an 
eigenenergy exactly equal to zero in the plot (the energy 
scale is referred to hwa + so that the plot shows only 
the non- linear part of the eigenspectrum). For ^ 
there is not such a possibility: the energies of the six 
dressed eigenstates are all non-zero (the fiatness of the 
dashed energy curve is just an effect of the large scale in 
the energy axis of the plot). The incoming photon can 
not leak into the cavity. Analogously, a photon of Ea can 
not penetrate into the cavity if a E"^^* photon is already 
there. 

The effects of the cavity dissipation and the atomic de- 
cay can be considered adding the term —ihT{a^a-\-wh) — 

i^lYl^j=i \ j) 0I> where F is the decay rate of the cavity 
and 7 is the atomic decay rate as in Section IIIII If we 
consider F,7 <C fid, however, their effects are to give a 
small linewidth to each eigenvalue, in any case insufficient 
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FIG. 5: (a) Eigenenergy for the doubly excited manifold of 
the system in Fig. 2J when the excitations are shared be- 
tween both the cavity modes (energies and gt in units of fid)- 
ga/^d = 0.3, y/fld = 0.01 and A/fld = O-l- The photon 
blockade effect with respect to a photon of colour uob is evi- 
dent if gb 7^ 0. (b) Restoration of photon blockade for a non- 
zero value of the single photon detuning 5 when two atoms 
are present in the cavity. The introduction of 5 drives the 
eigenstates of the doubly excited manifold out of resonance 
from the cavity field. No resonant transition is now found, 
showing the restored effectiveness of a blockade effect. The 
same parameters are used as in (a). 



to create suitable conditions for transitions to the dou- 
bly excited manifold. The external driving term induces 
transitions between the two levels in which the system 
has been effectively reduced. From the optical Bloch 
equations for a field interacting with a two-level atom 
(with a decay rate small compared to the Rabi frequency 
of the interaction) in stationary conditions, the popula- 
tion of the excited state of the atom is expected to reach 
^stat ^2/2 (the subscripts fc, I state the number of exci- 
tations in the ground and excited states). In our case, it is 
easy to prove that \£pump {down\ a \up)\ plays the role of 
a Rabi frequency for the interaction between the external 
field and the cavity-|-atom system, whose ground and ex- 
cited states are labeled as \down) and \up). For the cou- 
pling between \down) = |2,0,0) and \up) = \Da) 
with Enumv — O.irid (weak pump regime) and the pa- 
sta* _ 0.497. If the same 



calculation is performed with respect to the coupling be- 
tween \do'wn) = \Da) ® |0)^ and \Dab), which is the 
state corresponding to the dashed curve in Fig. (a) , we 
get crf2°* — 0.091 ^ 1/2. In the above discussion, we 
assumed the single-atom interaction. Let us consider a 
case in which a number of atoms are present in the cav- 
ity. This could happen because, experimentally, we do 
not have a perfect control of the number of atoms that 
cross the cavity. It has been shown in indeed that hav- 
ing more than a single atom in the cavity has the effect 
to introduce a large number of other energy levels that 
accumulate near the zero energy axis. This makes the 
photon blockade less effective. If we analyze our specific 
system, we see that the Hamiltonian itself acquires a col- 
lective operator structure. In general, atomic operators 
such as |i) (j| are replaced by Y^^i l*)fe '^i^h k label- 
ing the N particles inside the cavity. The injection of 
a photon into the cavity couples \'Q atoms ® \^) e ^^liZa) 
where we introduced the collective states: 



\D 



with 



ID 



atoms ®\^)e^ = 
\2) atoms 



|2-2), 



toml..atomN 



(g> |0) 



9a 11), 



mE„ (24) 



N 



{|122..2) + I212..2) 



I222..1)} 



atomL.atomN 



-^pump 

rameters in Fig. \^ we get 



being a singly excited symmetric Dicke state. To check 
the possibility for a photon blockade we have to seek 
the eigenvalues in the doubly excited manifold. Let us 
consider the simple case of iV = 2. The manifold with 
a single excitation is now composed of five bare states 
while they were three for the single atom case and there 
is an energy equal to huja [fujjb) if a photon from i?^^* 
{Ea^^) has excited the cavity field. Transitions to the 
singly excited manifold are possible. If we consider the 
doubly-excited subspace, we find eigenenergies suitable 
for a resonant transition that ruins the blockade effect. 
This behaviour can be verified for the cases of = 3, 4 
and we conjecture that this feature is present for an ar- 
bitrary TV. To bypass this problem, we introduce a single 
photon detuning in the transitions |2) ^ |4), |2) |6), 
while keeping the two-photon resonance. This implies the 
introduction of h5Y,j{\'^)j (4| + |6)j (6|) in the Hamilto- 
nian. If i5 > 7 these terms take account of the detuning 
of the dressed states off the cavity resonance: it shifts the 
energies in the doubly excited manifold and restores the 
blockade. The absence of resonant energies for N = 2, 
S and S ^ j is shown in Fig. O (b) . We have ex- 
tended this investigation up to the case of = 4. These 
results, obtained for values of the involved parameters 
achievable by current technology 25], are in agreement 
with the analysis in [2^ . 

We now introduce a quantum interference effect that 
explains the inhibition of transitions to highly excited 
states (the low values for cr^*'** can not be explained just 
by the detuning from resonance jH). For simplicity, we 
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consider the case of a single atom inside the cavity. If 
the external fields shine together the cavity, they simul- 
taneously try to excite the cavity field but, as long as one 
photon is fed into the cavity, there is no possibility for a 
second one to leak. We assume Ua < t^t and we include 
the two eigenstates, \Dab) and \D'^^), that in Fig. (a) 
are the nearest to the zero energy axis. This permits to 
estimate the blockade effect from a different perspective. 
The resulting effective five-level model is in Fig. |B1 Tak- 
ing a value of ga,gb ~ 0.557^, the splitting of \Dab) and 
\-^'ab) from the resonance becomes symmetric (as can be 
seen in Fig.[Sl (a), even for gb = 0.517^, ffa = 0.3V,d) and 
we take it to be A'. 

Pa'b) 



Probability 
1.0 




FIG. 6: Effective five-level model for the atom-f bi-modal 
cavity system, when the doubly excited eigenstates {\Dab) 
and \D'ab)) nearest to a resonant coupling with the single- 
quantum manifolds are included and we assume uja < uib- 
The detuning of \Dab), \D'ab) from LUa + u)b is symmetric for 
ga,gb = 5 — 0.5f2d. 

The effective Hamiltonian, in the basis composed of 
{|200) , \Da) , \Db) , \Dab) , \D'J} can be easily written for 
the couplings shown in Fig. 1^1 and its application to a 
state of the form 

\r,{t)) = Ci |200)+C2 \Da)+C^ \Db)+Ci \Dab)+C5 \D'J , 

(25) 

leads to a Schrodinger equation that can be recast in 
a set of five differential equations for the coefficients 
Cj (j = 1,..,5). We expect that the approximation 
of two-level system, consequence of the photon block- 
ade, has to be good, so that the populations of states 
\Dab) , l-D^b) should remain very small. Numerical inte- 
gration of these equations confirms our expectations and 
we get |C4|^ = IC'sl^ < 10^^. We thus neglect these terms 
and we find a simple analytical solution that results in 
Rabi oscillations of the populations of states |200), \Da), 

\Db) with the frequency flu = fld^pump/ + 5^ ' ^^^^ 
matches the analogous parameter from . The results 
are shown in Figs. [7| The amplitudes of oscillation for 
\C2\ and ICaT 

are 1/2, and they oscillate in phase (ab- 
solutely indistinguishable in the plot). This is because 
we do not know which one of them could have taken the 
excitation from outside. The probability of excitation is 
thus 1/2 for both the fields. At this time, the state of 
the system in an equally weighted superposition of \Da) 




Probability 
1.0 



FIG. 7: (a) Plot of the probabilities ICif (solid hne), \C2\ 
and ICal^ (dashed line). The latter two are completely su- 
perimposed, r is a rescaled time: r — t/^ln, with the 
effective Rabi frequency of the oscillations, (b) Comparison 
between the analytical (solid line: the same as in (a)) and the 
numerical (dashed line) solution for |Ci|^. The small offset 
of the numerical solution has been added to make the curve 
visible: the matching with the analytic solution is perfect. 



and \Di,). The calculation of C2 and C3 leads to: 
1 



V2 



{cos^J10>„, + cosz?b|01)„J® |2)^ 



1 



V2 



{sin^?a \l)atom + |3),,_} <E> |00),, 



(26) 



with cos I?, 



a, 6 



If we measure the state 



of the atom, as it exits from the cavity, and we find |2), 
the cavity modes are projected onto an entangled state 
with adjustable coefficients. Even if we have just roughly 
outlined the problem, this certainly deserves a further 
analysis. This result seems to be promising in the per- 
spective of quantum state engineering, in particular be- 
cause no single-atom addressing is required. 



VII. CONCLUSIONS 

The main result of this work is the discussion of the 
quantized picture of a model for double-EIT 13]. Our 
approach is based on a full Hamiltonian method that sim- 
plifies the problem of a many-level atomic system inter- 
acting with some electromagnetic fields. The results ob- 
tained have been confirmed by a dressed state approach. 
Here, we have shown the cross-phase modulation of two 
interacting fields. We suggested the Pr:YSO crystal as 
a candidate to embody the model. The success of the 
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quantization step has led us to investigate two particular 
problems in the context of quantum state engineering: 
the generation of entangled coherent states of light and 
a bi-chromatic photon blockade in CQED. In the former 
case, a scheme for the inference of the non-classicality of 
the generated state has been briefly discussed, including 
the effects of losses by the detection apparatus. 

For the photon blockade, we discussed the main fea- 
tures for the control of the population of two cavity field 
modes. Considering a flux of atoms crossing a bi-modal 
cavity, we have shown the effectiveness of the blockade 
effect even when more than a single atom is present in the 
cavity. With respect to a soHd state system (as Pr:YSO) 
placed into the cavity, the examined set-up is more suit- 
able for the realization of photon blockade because it 
overcomes a series of severe restrictions that in a multi- 
atom system are imposed because of the high-dispersion 



limit Q. 

In conclusion we have outlined specific applications of a 
particular kind of large and efficient non-linearity, namely 
the double-EIT regime, in the context of controlling a 
quantum system. 
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